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Summary 

A numerical model for the three dimensional thermal analysis of rocket thrust 
chambers and nozzles has been developed. The input to the model consists of the 
composition of the fuel/oxidant mixture and flow rates, chamber pressure, coolant 
entrance temperature and pressure, dimensions of the engine, materials and the 
number of nodes in different parts of the engine. The model allows for tempera- 
ture variation in three dimensions: axial, radial and circumferential directions and 
by implementing an iterative scheme, it provides nodal temperature distribution, 
rates of heat transfer, hot gas and coolant thermal and transport properties. 
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Nomenclature 


A 

C 

C p 

CCH 

ccw 

d 

D 

DCIN 

DG 

f 

FO% 

h 

i 

k 

L 

m 

P 

Pr 

r 

Re 

T 

THKNS 

V 

W 


area of the nozzle 

correction factor for heat transfer coefficient 

specific heat 

cooling channel height 

cooling channel width 

diameter of the thrust chamber 

characteristic diameter 

distance between bottom of the cooling 

channel and inner surface (see Fig. 3) 

gas side diameter 

friction factor 

fuel/oxidant percentage 

heat transfer coefficient 

enthalpy 

conductivity 

characteristics length 

total number of stations 

pressure 

Prandtl number 

radius 

Reynolds number 
temperature 
wall thickness 
velocity 
weight flow 
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X station position in longitudinal direction 

(distance from the rocket throat) 

Greek Symbols 

A S length of cooling channel between two stations 
p dynamic viscosity 

p density 


Subscrii 


A 

adiabatic 

C 

coolant 

G 

gas 

n 

station n 

s 

static 

w 

wall 

X 

reference 

0 

stagnation 
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Introduction 


The thermal analysis of spacecraft engines is an essential step in determining 
the rate of heat transfer from combustion gases to the coolant and in evaluating 
thermal stresses due to temperature gradients. In the past, rocket engines were 
designed to be used only once; therefore fatigue due to thermal stresses was not 
an issue arid there was no need for an accurate evaluation of temperature gra- 
dients. However, in the new reusable spacecraft engines such as SSME (Space 
Shuttle Main Engine), OTVE (Orbit Transfer Vehicle Engine) and IILLV (Heavy 
Lift. Launch Vehicle), reducing fatigue due to thermal stresses is an important fac- 
tor in increasing the life of the engine. 

The existing thermal model at NASA Lewis Research Center uses a uniform 
wall temperature assumption m both radial and circumferential directions. This 
assumption is reasonable as long as the hot gas pressure and temperature gradient 
are relatively small. However, in the new engines such as orbit transfer vehicle 
engines, the gas pressure and temperature gradients are relatively large and the 
chamber wall includes ribs, fins, and high LjD coolant passages, which make a uni- 
form wall temperature an unreasonable assumption. Existing thermal computer 
codes such as CINDA/SINDA [1], ADINAT [2], MARC [3], ANSYS [4]. and NAS- 
TRAN [5] cannot be used without major modifications because of the complex 
nature of the problem. The complexity of the problem is due to three-dimensional 
geometry, coolant, and hot gas heat, transfer coefficient dependence on the pres- 
sure and wall temperature, and unknown coolant pressure drop and properties. 
Therefore, a new thermal model, accounting for the complexities of the problem, 
needs to be developed. 
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The objective of this paper is to explain the numerical model developed to 
analyze rocket engine heat transfer characteristics. This model can be used to 
determine the temperature distribution in a regeneratively cooled thrust chamber 
by allowing for temperature variations in the radial, circumferential, and axial 
directions. 


Numerical Model 

The geometry of a typical regeneratively cooled thrust chamber is shown in 
Figure 1. The combustion chamber and nozzle wall consists of three layers made 
of ceramic coating, copper, and nickel with a number of cooling channels (layers 

of materials may change for different engines). The conductivities of the wall 
materials can be functions of the wall temperature. A numerical procedure for 
solving this problem is developed. In this procedure, the rocket thrust chamber 
and nozzle are subdivided into a number of two-dimensional stations along the lon- 
gitudinal direction (see Figure 2). The thermodynamic and transport properties 
of the combustion gases are evaluated using the computer program developed by 
Gordon and McBride [6-7]. The GASP [8] or WASP [9] programs are implemented 
to obtain coolant thermodynamics and transport properties. First the static pres- 
sures, temperatures and enthalpies for the combustion gases are evaluated using 
the ROCKET subroutine given in [6]. At a given station, by making an initial 
guess for the temperature distribution, the gas and coolant adiabatic wall temper- 
atures and wall properties can be evaluated based on the assumed wall temperature 
distribution using the computer properties codes [6-9] for combustion gases and 
coolant. The reference enthalpy of the gas side, tax is given by 

tGX = 0.5(?gvv + *gs) + 0.180(«go — *gs) (1) 
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where icw is a function of gas static pressure Pqs and gas side wall temperature 
Tgw and is evaluated using the program given in [6]. The gas side adiabatic wall 
enthalpy, Igaw is calculated using the following equation [10-11] 

IGAW = IGS + {PrGx) l * Z {iGO — * gs ) ( 2 ) 


where the gas reference Prandtl number Prcx is 


Ptgx — 


' Pax 


PGX 


&G X 


C Pax , hgx and kox are functions of Pgs and icx • The hot gas side heat transfer 
coefficient, he is given by 


ha = 

d G n 


( 3 ) 


where Cg„ is the gas side correlation coefficient and the Reynolds number is 
given by 


Rcgx = 


4 Wq Tgs 
ndG n fiGX Tgx 

Tgx = /{Pgs^gx) 

Tgs = f(PGS,*Gs) 

For the cooling channel the coolant velocity is calculated from the following 
equation: 


Vcs n = 


W c 


( 4 ) 


Pcs n Ac n N n 

Note that the coolant static density for the first iteration is set equal to the coolant 
static density of the previous station and at the first station it is set equal to its 
stagnation density. The static enthalpy is then calculated using the following 
equation: 

«cs„ = ico TT 2- ( 5 ) 
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The coolant stagnation enthalpy given in the right hand side of equation (5) is a 
function of the coolant stagnation temperature and pressure, i.e., 

ico! = }{Tco x i PcOi) (6) 


for the first station and for the other stations it is given by 

(<7n ~h Qn — 1 )^^n— l,n 


iCO„ = tCOn-i + 


2W C 


(7) 


where A5 n -i, n is the distance between two neighboring stations n — 1 and n which 
is calculated from 


A5 n _i, n = \/( ^ n Z 0 - 1 ) 2 + (*» - *n-l) 2 


( 3 ) 


and q n is the heat transferred from the hot gases to the coolant at station n. For the 
first iteration at station n, q n in equation (7) is not known; therefore the following 
equation is used to evaluate the stagnation enthalpy 


ico n = ico n -i + 


Qn — 1 ASn— l,w 

W c 


(9) 


The coolant static and reference Reynolds numbers respectively are given by 


•Rees = 


W c d Cn 

Ac n NnllCS 


( 10 ) 


and 

Recx = (11) 

Pcs Pew 

where ncs is a function of Pcs and ics and is calculated using the GASP program 
[8], or the WASP program [9] if the coolant is water. To employ a better value for 
the Reynolds number an average Reynolds number between the entrance and exit 
to each station is evaluated, i.e., 


R e cs Avg = 0.5(i?ecs„ + -RecSn-J 
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R e cx Avg = 0.5 {Recx n + 7?ecx n _, ) 

The friction factor in the cooling channel is calculated using the Moody diagram 
correlations assuming the cooling channel surface is smooth. The correlations for 
the friction factors are given by 


/ = 


— 64 


Recx 


/ = 4C, 


Avg. 


If Re CXAvg < 2.2 x 10 3 


.0014 + 


■125 




If Recx Avg > 2.2 x 10 3 & 


/ = .0778C a (Re C x Avg .)- 1021 


Recx Avg < 10 4 
If Recx Avg > 10 4 


( 12 ) 


where 


Ci = 


Cc 

.023 


= f (station) 


Once the friction factors are determined, the viscous pressure drop between stations 
n — 1 and n is calculated using Darcy’s law which is given by 

/ ( PCS„ + Pcs n _! 


i tergf ) + 

and the momentum pressure drop is calculated using the following equation: 


(13) 


( APcs„_i,„ )a/ 


' 2 1 

.(NA C )n-l+(VAc)n} 9c 


•[ 


{PCsAcN)n ( PCS^c 


1 ' 


(14) 


An average value of variables between stations n and n — 1 are used to improve the 
accuracy. The static pressure at each station is calculated based on the viscous 
and momentum pressure drops and is given by 


PcSn = PcSn., ~ [(APcs„_,, n )/ + ( APc5 n _! M (15) 

The reference and adiabatic wall enthalpies at the station are, respectively, calcu- 
lated from the following equations [10] 

icx n = 0.5(ics„ + icw n ) + 0.194(ico„ - ics n ) (10) 
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and 

icAWn — ics„ + (Prcx) l / 3 (ico„ - *cs„) (17) 

The adiabatic wall temperature is a function of the coolant static pressure and the 
adiabatic wall temperature and is evaluated using the GASP program [8]. Note 
that the Prandtl number in equation (17) is expressed by 

Pro X = 

kcx 

where C Vcx , pcx, kcx — f(Pcs,icx)- 


To evaluate the coolant heat transfer coefficient the correlation given in refer- 
ences [13-15] is used. In this correlation the Nusselt number is given by 

Nu 


Nu r 


= 0.023 Re 08 Pr 0A 


(18) 


where 

Nu r = V’ -0 ' 55 


ip = 1 + 0(Tw — T s ) 

and 

o \}_dp_\ _ 1 (fr)/ > 

^ l ~pdT ] ~ P ( JE )r 

Properties for the above correlation are based on the coolant static temperature 
Tcs, and static pressure Pcs ■ These properties are calculated using the GASP 
program [8], or the WASP program [9] if the coolant is water. It should also be 
noted that there are three coolant heat transfer coefficients and adiabatic wall 
temperatures. They are for the top, side, and bottom walls of the cooling channel. 
The variable heat transfer coefficient is due to the variable wall temperature in 
the cooling channel. The coolant reference and adiabatic wall enthalpies are also 
functions of wall temperature and are larger for the surface nodes closer to the 
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bottom of the cooling channel. 


Once the heat transfer coefficients and adiabatic wall temperatures for the 
hot gas and coolant are evaluated, the two-dimensional finite difference program 
discussed in reference [12] 1 2 3 is used to re-evaluate the wall temperature distribu- 
tion. Based on the revised wall temperature, new hot gas and coolant wall proper- 
ties, heat transfer coefficients and adiabatic wall temperatures are calculated using 
equations (1) through (18). Again, a new wall temperature distribution based on 
the most recent heat transfer coefficients and adiabatic wall temperatures is cal- 
culated using the two dimensional finite difference program [12]. This procedure 
is repeated until the relative difference between the temperature distribution of 
two consecutive iterations becomes negligibly small. By performing similar calcu- 
lations for all stations, the temperature profile for the entire engine is calculated. 
A complete flow chart of this numerical model is given in the Appendix. 

Three assumptions are made in the present model: 


1. Axial heat conduction through the engine wall is neglected, i.e., there is no 
heat conduction between two stations. 

2. Radiation heat transfer between the hot gases and surfaces of the engine is 
not included in the present analysis. 

3. Steady state flow and heat transfer. 

'The computer program discussed in [12] has been specifically developed for a rocket thrust 
chamber and nozzle two-dimensional configuration as shown in Figure 1 . Because of the symmetry 
of the configuration, computations are performed for only one cell (see Figure 3). In this program 
the thermal conductivities can be functions of temperature and it implements a successive over- 
relaxation formula for a quick convergence. The resuling output of the program consists of the 
nodal temperature distribution, heat transfer to the coolant and heat transfer from the hot gas. 


10 




A discussion regarding the validity of the above assumptions is presented in 
the next section. 


Results and Discussion 


A computer code based on the numerical model introduced in the previous sec- 
tion has been developed. A complete description of the code is presented in [16] 
and it is used to determine the temperature distribution and heat transfer char- 
acteristics of a Liquid Oxygen/Kerosene rocket engine (engine 700). The engine 
has the following specifications: 


Fuel 

Oxident 

Coolant 

Chamber stagnation pressure Pqo 
Coolant stagnation pressure Pco 
Fuel weight flow rate 
Coolant weight flow rate 
Fuel/Oxident percentage 
Coolant stagnation temperature 
Number of cooling channels 


Kerosene (RP-1) 
L0 2 

lo 2 

2000 psi 
5000 psi 
58.023 lb/sec 
39.2 lb/sec 
26.316% 

150° R 
100 


This engine is subdivided into 22 stations. Table 1 shows dimensions of the 
engine and some thermal characteristics at each station (see Figure 3 for notation). 
Note that dimensions given in Table 1 are in inches. There are two variables that 
remain constant along the engine. These varables are: Cc = 5.36 x 10~ 4 and 
DCIN = 0.035 in. The outer surface radiates to space and its emissivity is 0.9. 
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The thermal conductivities of wall materials, i.e., nickel and copper are functions 
of temperature and are given by Figures 4 and 5 [19], and no coating is used in the 
engine. Tables 2, 3 and 4 show the resulting temperature distribution for stations 
1, 9(throat) and 22 respectively. The previous model at NASA Lewis Research 
Center used a lumped temperature assumption in the radial and circumferential 
directions while the present model allows temperature variations in all three di- 
rections. Close examination of the temperature distributions given in Tables 2, 3 
and 4 reveals that the temperature gradient is relatively large in radial direction, 
especially for station 9. Hence, the uniform wall temperature assumption in the 
radial and circumferential direction used in the previous NASA model is not a 
good assumption for the engine considered here. This may also be true for any 
other high pressure thrust chamber. 

Figures 6 and 7 show temperatures of the inner and outer surfaces and the 
coolant along the engine for chamber pressures of 2000 and 1200 psi respectively. 
The resulting lumped wall (1-Dimensional model) temperatures along radial and 
circumferntial directions are also given in Figures 6 and 7. As can be seen from 
these figures, the lumped temperature results are very close to those of the inner 
wall at the early stations; however, in the later stations the difference becomes 
larger which might be due to the large temperature gradient. For engines with 
lower chamber pressures, the temperature gradient is small; hence, the lumped 
temperature assumption will provide reasonable results. 

Despite the improvements of the present code over the existing model, some 
items remain that should be considered in future modifications of the present 
model. One item is the axial heat conduction through the cooling channel wall. 
By looking at the wall temperature distribution along the longitudinal direction 
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given in Figures 6 and 7, one can see that there is a large wall temperature gradient 
in the axial direction. Hence, some heat is transferred by conduction between two 
neighboring stations. This heat transfer can roughly be estimated through 


9n— l,n — 


kA n — l,n 
A'S'n-l.n 


(Tn-1 


~ T n ) 


(19) 


where < 7 „_i in , A n - i iU and are the heat conduction, contact area and dis- 

tance between stations n — 1 and n respectively. 


Another item which needs serious consideration is the radiative heat transfer 
from the hot gas to the surface and from surface to surface. The chamber tempera- 
ture for the engine considered here reaches 3800° R while the wall temperature is of 
the order of 1200°R. With this temperature range, radiation becomes a significant 
means of heat transfer. A rough estimate of the radiative heat transfer from the 
hot gas to the surface can be obtained from [17-18] 

q G ^s = I<t( 1 - wo )GS°{T* - Tj) (20) 

where Tq and Ts are gas and surface temperatures respectively, K t and uq are 
gas total extinction coefficient and albedo for scattering, GS total exchange factor 
from gas to surface and a is the Stefan Boltzmann constant. Besides radiative 
transfer between the hot gas and surface, there is some radiative transfer from the 
high temperature surface to the lower temperature surface. 

Finally, most rocket engines perform for a relatively short period of time and, 
for maneuvering purposes, the fuel and coolant flow rate vary during this. There- 
fore, they never reach a steady state condition. A transient thermal analysis, 
although a bit ambitious, would make the present model capable of simulating 
actual rocket performance. 
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Concluding Remarks 


A numerical model for three dimensional thermal analysis of rocket engines has 
been developed. This model allows temperature variation along three directions: 
axial, radial and circumferential. The numerical results presented show that there 
is a large temperature gradient in the axial direction for engines with a high cham- 
ber pressure. This makes the lumped temperature assumption used in the previous 
model invalid for these types of engines. 

Experimental results exist for the engine analyzed in this paper. However, a 
comparison of the experimental and numerical results is not possible for two rea- 
sons. First, the engine only ran for five seconds during the experiment and did 
not reach a steady state condition. The fuel weight flow rate varied during this 
period. The second reason is that the exact position of the thermocouples used for 
the experiment is not known. A more carefully controlled experiment is needed to 
make a comparison with the numerical results meaningful. 

The numerical model needs to be modified further to incorporate a transient 
analysis, effects of axial conduction and chamber radiation. Efforts are presently 
under way to include these items into the numerical model. 
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n = 1 


n = n — 1 


MAKE AN INITIAL GUESS FOR THE WALL 
TEMPERATURE DISTRIBUTION AT STATION n 


CALCULATE REFERENCE 
ENTHALPY*, i GX , EQ. (1) 


CALCULATE ADIABATIC WALL 
ENTHALPY*, i GAW , EQ. (2) 


CALCULATE GAS ADIABATIC WALL 
TEMPERATURE* T GA w = f(Pas,ia A w) 


CALCULATE HOT GAS SIDE HEAT 
TRANSFER COEFFICIENT h G EQ. (18) 
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© 


n — 1 


n > 1 


Pcs i — pco 
*cOj = ico 


Pcs i = Pco 


PCSn = / , cs (n _ 1) 


j > 1 


Pcs„ — f (Pcs\ > i J cs n ) 


CALCULATE COOLANT STAGNATION ENTHALPY 
EQ. (7) WHEN j > 1 AND EQ. (9) WHEN j = 1 


CALCULATE COOLANT VELOCITY EQ. (4) 


CALCULATE STATIC ENTHALPY EQ. (5) 


CALCULATE REYNOLDS NUMBERS Re cs , Recx 
AND AVERAGE REYNOLDS NUMBERS 
EQ. (10) AND EQ. (11) 


CALCULATE FRICTION FACTOR / EQ. (12) 
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COMPUTE VISCOUS PRESSURE DROP EQ. (13) 


FIND MOMENTUM PRESSURE DROP EQ. (14) 


EVALUATE STATIC PRESSURE EQ. (15) 


FIND COOLANT WALL PROPERTIES* 
Cp cw ,Pcw, k C w, Pew, icw = f (Pcs, Tew) 


CALCULATE COOLANT REFERENCE 
ENTHALPY i cx EQ. (16) 


FIND COOLANT REFERENCE PROPERTIES+ 
Cp cx ,Pcx, kcx,Pcx,icx = / (Pcs, Tex) 


FIND COOLANT STATIC PROPERTIES* 
Cp cs ,fics,kcs,Pcs,ics = f(Pcs,Tcs) 


FIND COOLANT ADIABATIC ENTHALPY, 
icAW EQ. (17) AND ADIABATIC WALL 
TEMPERATURE* T C aw = f{Pcs,ic A w ) 
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"Use the computer program given in [61. 

+ Use the computer program given in [o] and [9]. 
$ Use the computer program given in [12]. 
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Table 1: Parameters of thrust chamber and nozzle at different stations. 


(Dimensions are in inches) 


Station 

X 

DG 

CCW 

CCH 

THKNS 

Cg 

1 

7. 

6.26 

0.06 

0.191 

0.462 

0.0234 

2 

5.024 

5.2 

0.06 

0.186 

0.423 

0.0212 

3 

3.5 

4.382 

0.06 

0.185 

0.42 

0.0185 

4 

2. 

3.58 

0.06 

0.187 

0.422 

0.0168 

5 

1 . 

3.042 

0.04 

0.16 

0.395 

0.0162 

6 

0.5 

2.778 

0.04 

0.146 

0.381 

0.0143 

7 

0.337 

2.688 

0.04 

0.135 

0.37 

0.0144 

8 

0.237 

2.642 

0.04 

0.124 

0.359 

0.0152 

9 

0 . 

2.6 

0.04 

0.095 

0.33 

0.0165 

10 

-0.274 

2.656 

0.04 

0.087 

0.322 

0.0168 

11 

-0.506 

2.7 

0.04 

0.077 

0.312 

0.0192 

12 

-0.906 

2.924 

0.04 

0.06 

0.295 

0.0219 

13 

-1.306 

3.092 

0.04 

0.066 

0.301 

0.022 

14 

-1.506 

3.178 

0.04 

0.069 

0.304 

0.022 

15 

-1.906 

3.344 

0.04 

0.071 

0.306 

0.0225 

16 

-2.906 

3.77 

0.04 

0.078 

0.313 

0.023 

17 

-3.506 

4.022 

0.06 

0.077 

0.312 

0.023 

18 

-4.906 

4.468 

0.06 

0.095 

0.33 

0.023 

19 

-5.906 

4.666 

0.06 

0.112 

0.347 

0.023 

20 

-6.906 

4.774 

0.06 

0.112 

0.347 

0.023 

21 

-7.572 

4.8 

0.06 

0.11 

0.345 

0.023 

22 

-10. 

4.8 

0.06 

0.11 

0.345 

0.023 
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Table 2: Temperature distribution ( °R ) at station 1. 



RADIANS FROM CENTERLINE OF CHANNEL 

| RADIUS (in) 

0.03141 

0.02410 

0.01679 

0.00948 

0.00632 

0.00316 

0.00000 


221.16 

221.12 

221.02 

220.91 

220.87 

220.85 

220.84 


221.31 

221.23 

221.04 

220.82 

220.75 

220.71 

220.69 

3.4560 

222.04 

221.80 

221.16 

220.39 

220.14 

219.98 

219.92 

3.4060 

224.72 

223.95 

221.81 

218.86 

217.76 

217.01 

216.75 

3.3560 

233.84 

231.86 

225.83 

215.20 

207.78 

203.79 

202.53 

3.3178 

243.40 

241.37 

235.37 

225.77 




3.2796 

263.25 

260.94 

254.06 

242.82 


COOLING 


3.2414 

295.26 

292.54 

284.25 

270.18 


CHANNEL 


3.2032 

342.21 

339.35 

330.28 

313.28 




3.1650 

405.61 

403.93 

398.86 

389.60 

402.16 

406.84 

408.12 

3.1562 

421.65 

420.25 

416.43 

412.42 

417.54 

420.72 

421.75 

3.1475 

438.18 

437.02 

434.13 

432.22 

434.49 

436.36 

437.03 

3.1388 

455.15 

454.14 

451.80 

450.69 

451.87 

453.01 

453.46 

3.1300 

472.52 

471.56 

469.39 

468.50 

469.39 

470.30 

470.67 
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Table 3: Temperature distribution at station 9(throat). 



RADIANS FROM CENTERLINE OF CHANNEL 

RADIUS (in) 

0.03141 

0.02410 

0.01679 

0.00948 

0.00632 

0.00316 

0.00000 

1.6300 

445.83 

445.83 

445.83 

445.83 

445.83 

445.82 

445.82 

1.5800 

445.84 

445.84 

445.83 

445.82 

445.82 

445.81 

445.81 

1.5300 

446.04 

446.01 

445.92 

445.81 

445.71 

445.64 

445.61 

1.4800 

448.38 

448.05 

447.12 

445.76 

444.40 

443.41 

443.06 

1.4300 

477.19 

474.38 

465.80 

450.96 

426.24 

412.99 

408.81 

1.4110 

504.65 

502.54 

496.36 

486.62 




1.3920 

558.36 

556.06 

549.23 

538.01 


COOLING 


1.3730 

641.91 

639.14 

630.78 

616.72 


CHANNEL 


1.3540 

763.19 

760.31 

751.22 

734.35 




1.3350 

926.10 

926.25 

927.01 

929.35 

973.30 

992.40 

997.89 

1.3262 

1002.06 

1002.96 

1006.17 

1013.56 

1028.32 

1038.54 

1042.02 

1.3175 

1076.43 

1077.36 

1080.31 

1085.55 

1092.04 

1097.06 

1098.89 

1.3087 

1149.03 

1149.81 

1152.09 

1155.60 

1159.20 

1161.98 

1163.01 

1.3000 

1220.29 

1220.99 

1222.99 

1225.92 

1228.76 

1230.91 

1231.71 
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Table 4: Temperature distribution at station 22. 



RADIANS FROM CENTERLINE OF CHANNEL 

RADIUS (in) 

0.03141 

0.02410 

0.01679 

0.00948 

0.00632 

0.00316 

0.00000 

2.7450 

869.30 

869.25 

869.13 

868.98 

868.89 

868.83 

868.81 

2.6950 

869.73 

869.61 

869.29 

868.89 

868.67 

868.52 

868.46 

2.6450 

872.20 

871.65 

870.17 

868.18 

867.02 

866.21 

865.91 

2.5950 

884.41 

882.13 

875.52 

865.38 

858.03 

852.75 

850.85 

2.5450 

939.58 

933.48 

914.10 

877.00 

803.89 

766.13 

754.57 

2.5230 

957.83 

952.86 

938.27 

915.86 




2.5010 

995.05 

990.63 

977.55 

956.39 


COOLING 


2.4790 

1049.73 

1045.61 

1032.95 

1010.88 


CHANNEL 


2.4570 

1121.36 

1118.15 

1107.72 

1086.76 




2.4350 

1207.01 

1206.33 

1204.84 

1204.56 

1236.81 

1249.62 

1253.25 

2.4262 

1242.25 

1242.30 

1243.28 

1248.95 

1263.71 

1272.92 

1275.90 

2.4175 

1277.70 

1278.22 

1280.50 

1286.94 

1295.09 

1301.16 

1303.31 

2.4087 

1313.06 

1313.82 

1316.64 

1322.78 

1328.31 

1332.65 

1334.25 

2.4000 

1348.17 

1349.00 

1351.96 

1357.89 

1362.71 

1366.50 

1367.92 
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FIGURE 1. - A ROCKET THRUST CHAMBER AND NOZZLE. 
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TEMPERATURE, °R 

FIGURE A. - THERMAL CONDUCTIVITY OF NICKEL. 



0 400 800 1200 1600 2000 2400 2800 3200 3600 

TEMPERATURE, °R 

FIGURE 5. - THERMAL CONDUCTIVITY OF COPPER. 
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FIGURE 6. - INNER AND OUTER WALL SURFACE AND COOLANT 
TEMPERATURE WHEN P G0 = 2000 PS I AND P CQ = 5000 PS I. 
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TEMPERATURE 



FIGURE 7. - INNER AND OUTER WALL SURFACE AND COOLANT 
TEMPERATURE WHEN Pg 0 = 1200 PSI AND P CQ = 3200 PSI. 
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